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Abstract. Using a sample of unprecedented size (about 400 objects) of H ll galaxies in which the oxygen abun- 
dances have been obtained using the temperature derived from the [O ill] A4363/5007 line ratio, we confirm that 
the H II galaxies form a very narrow sequence in many diagrams relating line ratios and H/3 equivalent width. We 
divide our sample in three metallicity bins, each of which is compared with sequences of photoionization models 
CO , for evolving starbursts with corresponding metallicity. Our aim is to find under what conditions a theoretical 

sequence can reproduce all the observed trends. Taking into account the presence of an older, non-ionizing stellar 
population, for which independent indications exist, we find that the simple model of an adiabatic expanding 
bubble reproduces the observational diagrams very well if account is taken of an aperture correction and the 
covering factor is assumed to decrease with time exponentially with an e-folding time of 3 Myr. We find that the 
He II A4686 nebular line emission occurs too frequently and in too wide a range of EW{Hl3) to be attributable 
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. to either the hard radiation field from Wolf-Rayet stars or the X-rays produced by the latest stellar generation. 



' Assuming that the He ll A4686 line is due to photoionization by a hot plasma at a temperature of 10® K, a 

^SJ , total X-ray luminosity of 10*° — 4 x 10*° erg s~^ is required for at least half of the sources. We find evidence for 

■ self-enrichment in nitrogen on a time scale of several Myr, and argue for a possible self-enrichment in oxygen as 

well. 
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^ 1. Introduction 1989; McGaugh 1991; Cid Fernandes et al. 1992; Bernlohr 

1993; Cervino & Mas-Hesse 1994; Garcia- Vargas & Diaz 

^> It has been known for a long time that giant H ii regions ^994. Garcia- Vargas et al. 1995; Stasinska & Leitherer 

k> form a well defined sequence m emission hne diagnostic di- ^995 (hereafter SL96); Bresolin et al. 1999; Dopita et al. 

^ agrams (McCall et al. 1985; VeiUeux & Osterbrock 1987). 20OO; Moy et al. 2001; Chariot & Longhetti 2001; Stasinska 

, The first interpretations of this sequence were based on 2OOI (hereafter SSLOl); Zackrisson et al. 2001). 

single star photoionization models, and it was concluded ^hese studies used different prescriptions for the modeling 

that the driving parameter of the sequence is the metallic- ^f the stellar population as well as of the nebular emission, 

ity, and that variations of an additional parameter - either Comparisons with observations are difficult because there 

the effective temperature of the ionizing stars or the ion- j^^^^ independent parameters that drive the 

ization parameter ~ are required in order to reproduce the emission line properties of extragalactic H 11 regions: age 

observed sequence (McCall et al. 1985; Dopita & Evans metallicity. The general conclusion from these stud- 

^^^^)- ies is that the Salpeter initial mass function with an up- 

However, giant H il regions are powered by intensive per stellar mass limit around 100 M© is consistent with 

bursts of star formation (e.g. Sargent & Searle 1970; Mas- the observed emission line ratios (although Bresohn et al. 

Hesse & Kunth 1991). Stellar evolution produces a gradual (1999) argue for a lower cut-off mass at high metallicity). 

change with time of the integrated stellar energy distribu- However, several problems are noted. All the models for 

tion, which has to be taken into account in the modeling instantaneous starbursts predict a strong drop in [O III] 

of giant H 11 regions. Numerous studies have produced A5007/H/3 at about 5 Myr. This corresponds to the max- 

grids of photoionization models for H 11 regions consider- imum lifetime of O stars so the prediction is very robust, 

ing that aspect (e.g. Tcrlcvich & Mclnick 1985; Olofsson The age of a starburst in an isolated giant H II region can 

be estimated, to a first approximation, from the equivalent 

Send offprint requests to: G. Stasiriska, width of H/3, £'IF(H/9) (this is not feasible in giant H II re- 
e-mail: grazyna. stasinska@obspm.fr 
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gions belonging to spiral galaxies, where old stellar popu- 
lations strongly contribute to the continuum at H/3). Using 
samples of isolated extragalactic H ii regions, hereinafter 
referred to as H ii galaxies, SSLOl and Zackrisson et al. 
(2001) showed that the observed drop was rather mild 
and displaced towards lower values of EW{llf3). SSLOl 
attributed this to the effect of underlying old populations 
of stars, shown by Raimann et al. (2000) to be present 
in H II galaxies. Other causes such as leakage of ioniz- 
ing photons or selective dust extinction were also men- 
tioned. The classical diagnostic diagram [O ill] A5007/H/3 
vs. [O ii] A3727/H/3 does not seem to be completely under- 
stood in terms of pure photoionization models. One way 
out is to advocate some additional heating (SSLOl) or the 
contribution of zones of low ionization parameter (Moy 
et al. 2001). A very significant trend of [O i] A6300/H/3 
increasing with decreasing EW(lif3), discovered by SL96 
and SSLOl remains to be quantitatively explained. 

The purpose of the present paper is to quantify the 
conditions needed to reproduce the observed emission line 
sequence of giant H II regions, including the very small 
dispersion seen in many emission line diagrams. For a more 
meaningful comparison of models with observed data, our 
observational sample is composed of objects for which the 
age of the ionizing star cluster and the metallicity can be 
estimated in an independent way. This limits the sample to 
H II galaxies with measured [O ill] A4363 line intensities. It 
is only when such a sample is fully understood that one has 
a chance to better understand the entire extragalactic H ii 
region sequence, including more metal-rich H li regions 
such as those found in the inner parts of disk galaxies. 

In Sect. 2 we describe the observational sample to 
which we compare the models and the preliminary treat- 
ments of the data. In Sect. 3 we outline the modeling strat- 
egy, and in Sect. 4 we present our results. The main find- 
ings of our paper are summarized and discussed in Sect. 
5. 



2. The observational sample 

The total sample consists of two subsamples. 

Galaxies in the first subsample have been extracted 
mainly from the First and Second Byurakan surveys. The 
total number of the galaxies is ~ 100. The high signal- 
to-noise ratio spectra of these galaxies in the wavelength 
range A3600 - 7400A have been obtained with different 
2m - 10m class telescopes. Some galaxies were observed 
several times. We included in the sample only one obser- 
vation of each galaxy for which the spectrum was obtained 
with the highest signal-to-noise ratio. The spectra were 
reduced in a homogeneous way according to prescriptions 
by Izotov et al. (1994, 1997) and have been used for He 
and heavy element abundance determination in a series 
of published papers. The line intensities corrected for ex- 
tinction can be found in Izotov et al. (1994, 1997, 1999, 
2001a, 2001b), Izotov & Thuan (1998), Thuan et al. (1995, 
1999), Lipovetsky et al. (1999), Guseva et al. (2000, 2001). 



Hereafter we will call this subsample the "Byurakan sam- 
ple" . 

The second subsample includes ~ 300 emission-line 
galaxies from the early data release (EDR) of the 
SLOAN digital sky survey (SDSS) (Stoughton et al. 2002). 
Hereafter we will call this subsample the "SDSS sample" . 
The flux-calibrated spectra of these galaxies in the wave- 
length range ~ A3820 - 9300A obtained with 3" round 
slits have been extracted from the Space Telescope Science 
Institute archivesl^. As in the case of the Byurakan sam- 
ple we included in the SDSS sample only galaxies with a 
detected [O iii] A4363 emission line. The extracted spec- 
tra have been transformed to the linear wavelength scale 
and zero redshift. The line intensities have been measured 
by Gaussian fitting to the line profiles. The errors in line 
intensities include the errors of Gaussian fitting and those 
of the continuum level. Then the line intensities have been 
corrected for the interstellar extinction and underlying hy- 
drogen stellar absorption lines from the observed Balmer 
decrement. The emission line intensities relative to H/3 
and the H/3 equivalent widths are available on request to 
Y. Izotov. Because the spectra in the SDSS sample start 
from ^ A3820A the hue [O ii] A3727 has not been detected 
in low-redshift galaxies with z ^ 0.025. In those cases 
the intensity of the [O ii] A3727 emission line was calcu- 
lated from the intensities of the [O ii] A7320,7330 emission 
lines. Another problem with the SDSS spectra is that in 
many cases some strongest lines, most often Ha and [O ill] 
A5007, are clipped. Therefore, many objects with strong 
emission lines were not included in our sample. However, 
in the cases when only Ha and [O ill] A5007 emission lines 
are clipped, their intensities are adopted to be respectively 
2.8x/(H/3) and 3x/([0 iii] A4959). In general, the galax- 
ies from the SDSS sample are fainter than those from the 
Byurakan sample and their spectra are of lower signal-to- 
noise ratio. 

The main characteristics of the Byurakan and SDSS 
samples are compared in Fig. 1 which shows histograms 
of the values of redshifts z (panel a), H/3 equivalent widths 
EW{Ilf3) (panel b), colour excesses E{B - V) as de- 
rived from the Balmer decrement (panel c) and H/3 lu- 
minosities, i(H/3) (panel d). SDSS sample data are rep- 
resented by a solid line and Byurakan sample data by 
a dotted line. On average, the SDSS galaxies are more 
distant but the H/3 luminosities are similar to those of 
the Byurakan sample. The equivalent width distributions 
are different, being more skewed toward large equivalent 
widths in the Byurakan sample. Note that many galaxies 
from the Byurakan sample have been selected for the pri- 
mordial helium abundance determination, favouring large 
equivalent widths of the emission lines. On the other hand, 
since the SDSS galaxies are more distant on average, a 
larger fraction of galaxy light falls within the slit, con- 
tributing to the continuum at H/3. However, wc found no 
correlation between EW(R(3) and redshift, which seems 
to indicate that statistically, the light from the underly- 

^ http://archive.stsci.edu/sdss. 
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ing galaxy does not contribute substantially to EW{Hf3). 
The distributions of E{B — V) are roughly similar in both 
samples. For the majority of the objects, the extinction is 
small. 



2.1. Plasma diagnostics 

The electron number densities derived from [S ii] 
A6731/6717 line ratios in the objects from the Byurakan 
sample are small (around 100 cm"'^ or smaller for most 
objects). They are more uncertain in the galaxies from 
the SDSS sample because of the lower signal-to-noise ra- 
tio spectra. Therefore, in the abundance determinations, 
a density of 100 cm"'^ was assumed for those objects. The 
oxygen abundances have been obtained from the redden- 
ing corrected emission line intensities with the electron- 
temperature-based method. The temperature T[o iii] de- 
rived from [O III] A4363/5007 was taken as representa- 
tive of the [O III] A5007 and H/? emission, and the [O ii] 
A3727 line was assumed to be emitted at a temperature of 
0.7x(r[o III] - 10,000 K) + 10,000 K as in SL96. The ion- 
ization correction for ionization stages higher than O'*'^ 
is minimal. The derived oxygen abundances in our sample 
range between 10^^ and 2 x lO^'*. 

2.2. Comparison to former samples 

Our merged sample is far more numerous than the Izotov 
or Popescu & Hopp (2000) samples used in SSLOl, since 
it contains around 400 objects. It is of comparable size 
to the entire sample of H ii galaxies from the Terlevich 
et al. (1991) survey, which comprised about 300 objects. 
However, being limited to objects with low metallicities 
and EW{H0) larger than lOA due to the condition im- 
posed on the detection of [O ill] A4363, our sample is four 
times more numerous than the Terlevich sample per bin of 
metallicity and equivalent width. In addition, the signal- 
to-noise ratio in our sample is generally much higher than 
in the Terlevich sample. This gives us the opportunity to 
investigate the evolution of H ii galaxies using statistical 
considerations with unprecedented detail. 

2.3. Observational trends 

Figure 2 presents our entire sample of H ii galaxies in a set 
of diagrams involving emission line ratios and EW(Hl3). 
The number of data points present in each diagram is 
shown in the upper right of each panel. The trends noted 
in SSLOl are clearly confirmed. Most conspicuous are the 
decrease of [O iii] A5007/[O ii] A3727 and the increase of 
[O i] A6300/([O II] X3727+ [O iii] A5007) as EW{lll3) 
decreases, as well as the constancy of He i A5876/H/3 
and ([O ii] A3727 + [O iii] A5007)/H/3 and the mild de- 
crease of [O III] A5007/H/3. There is also a clear tendency 
for [N ii] A6584/[0 ii] A3727 to increase as EW(Jil3) de- 
creases, although the relation is more dispersed. The pure 
emission line diagnostic diagrams [O ill] A5007/H/? vs. 



[O II] A3727/H^ and [O ill] A5007/H/? vs. [O i] A6300/H/? 
do not draw such an extended sequence as when using 
data including objects with higher metallicities, such as 
H II regions located in the central parts of spiral galaxies 
(McCaU et al. 1985; van Zee et al. 1998a), but the points 
are strongly gathered in a well defined zone of the plane. 
When plotting the Byurakan sample and the SDSS sam- 
ple separately, the points are indistinguishable in these 
diagrams in the zone of overlap of EW{ll(3), except that 
the He i A5876/H/3 shows much less scatter due to the 
higher quality of the data. This indicates that the nature 
of the objects is similar and justifies the merging of the 
two samples for the study of the emission line sequence in 
H II galaxies. 

In the He ii A4686/H/3 vs. EW(Rl3) diagram, the num- 
ber of points is much smaller than in the previous dia- 
grams. The reason is not only that He ii A4686 is a weak 
line, but also that this line is difficult to measure since it is 
superimposed on the Wolf-Rayet bump seen in the spectra 
of many H ii galaxies (Guseva et al. 2000). It could not 
be extracted in the fainter objects from the SDSS sample 
because of the lower signal-to-noise ratio in their spectra, 
and could be safely measured only in 48 objects from the 
Byurakan sample. The lower envelope is obviously an ob- 
servational effect. What is important to note is that the 
dispersion is rather large, that nebular He II A4686 line 
is detected in the entire range of EW{}if3) of our sam- 
ple, and that no tendency is seen with EW(ll(3) (see also 
Guseva et al. 2000). 

The fact that a sample of 400 objects shows such a con- 
spicuous collective behaviour in the plots of Fig. 2 is re- 
markable, given that detailed studies of a few of them show 
that their nature is rather complex (Thuan et al. 1997, 
1999; Papaderos et al. 1998, 1999; Stasiriska & Schaerer 
1999). This means that, in spite of this complexity, the 
evolution of H ii galaxies obeys simple laws. In the follow- 
ing, we will aim at describing this evolution in an empirical 
way using photoionization models, in order to find some 
clues on the driving mechanisms of this evolution. 

3. Modeling 

3.1. The modeling codes 

The computations are done with the same codes as used 
by SSLOl. The theoretical spectral energy distributions 
are obtained from the evolutionary synthesis models of 
Schaerer & Vacca (1998)0. They use the non-rotating 
Geneva stellar evolution models, with the high mass- loss 
tracks of Meynet et al. (1994). The spectral energy dis- 
tributions for massive main-sequence stars are taken from 
the CoStar models (Schaerer & de Koter 1997), which in- 
clude the effects of stellar winds, non-LTE, and line blan- 
keting in an approximate manner. The pure He models of 
Schmutz et al. (1992) are used for Wolf-Rayet stars. The 
spectral energy distributions from the plane-parallel LTE 

^ The predic tions from these synthesis models are available 
on the Web at http://webast.ast.obs-mip.fr/people/schaerer/. 
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Fig. 1. Histograms of rcdshifts z (panel a), H/3 equivalent widths (panel b), colour excesses E(B — V) as derived from 
the Balmer decrement (panel c) and H/? luminosities (panel d). SDSS sample data are represented by a solid line and 
Byurakan sample data by a dotted line. 



models of Kurucz (1991) are used for the low mass stars 
which build up the continuum. 

Recent computations of stellar energy distributions us- 
ing more complete stellar atmosphere codes (Hillicr & 
Miller 1998; Pauldrach et al. 2001) show significant dif- 
ferences with the models used in the present work. The 
most important differences, however, concern high mctal- 
licity stars. At the metallicities relevant for our studies, 
the agreement with the atmospheres used here is reason- 
able. In general, the ionizing radiation predicted by the 
new model atmospheres is softer than that predicted by 
the model atmospheres entering the population synthesis 
of Schaerer & Vacca (1998). This concerns both the main 
sequence 0-stars (resulting in a smaller ratio of He° to H° 
ionizing photons QHe/Qn) and Wolf-Rayet stars (result- 
ing in a smaller Qhc+ZQu)- Stellar population synthesis 
models including these new atmosphere models were not 
yet available at the time when the present work was done. 
But they are not expected to significantly alter the con- 
clusions of the present work. 

The photoionization models using the spectral energy 
distributions as an input are built with the code PHOTO, 
as described in SL96. 



3.2. The modeling strategy 

Instead of producing wide grids of models, as in the stud- 
ies listed in the introduction, we try to fit the observed 



data with sequences of models that would represent the 
true evolution of a giant H ii region excited by a cluster 
of evolving stars. We take advantage of the fact that we 
have determined the oxygen abundance (hereafter sim- 
ply referred to as the metallicity) for all the objects of 
our sample, so that each object can be compared to one 
model of given metallicity. To make the problem visually 
tractable, we divide our sample of objects in three metal- 
licity bins: the "high" metallicity bin contains all the ob- 
jects of our sample having 0/H > lO"** ( i.e. larger than 
0.12 times the solar oxygen abundance as given by Anders 
& Grevesse 1989), the "intermediate" one contains all the 
objects with 0/H between 10~^ and 3x 10~^ ( i.e. between 
0.12 and 0.036 times solar) and the "low" metallicity bin 
contains all the objects with 0/H < 3 x 10~^ ( i.e. less 
than 0.036 times solar). Obviously, the words "high", "in- 
termediate" and "low" are not to be taken literally in this 
context, since our entire sample is composed of objects 
with 0/H < 2 X lO""* (i.e. less than 0.25 times the solar 
value). We have not included in any metallicity bin the ob- 
jects (from the SDSS data) for which the estimated error 
in the [O ill] A4363 intensity is larger than 50% in order 
to limit undue shift of metallicity bins as a consequence of 
important errors in the abundance determinations. Each 
metallicity bin is compared to models with corresponding 
metallicity: Z = 0.2 Z©, Z = 0.05 Zq, and Z = 0.02 Zq, 
for the "high" , "intermediate" and "low" metallicity bins 
respectively. As in SSLOl and in most recent studies on 
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Fig. 2. The entire sample of H II galaxies in a set of diagrams involving emission line ratios and EW(Rl3). The 
number of data points present in each diagram is shown in the upper right of each panel. 



extragalactic H ii regions, we adopt as a reference the solar 
abundances from Anders & Grevesse (1989), although the 
latest determinations of solar abundances are significantly 
different. Especially, the oxygen abundance in the Sun as 
derived by AUende Prieto et al. (2001) is 4.9 x 10~* instead 
of 8.5 X 10~"^ by Anders & Grevesse (1989). In the present 
work, the term "solar" is used for an easier reference with 
previous work. As in SSLOl, we use the McGaugh (1991) 
prescription to link the abundances of all the elements to 
the abundance of oxygen. In particular, He/H = 0.0772 + 
15 (0/H) and log N/0 = 0.5 log (0/H) + 0.4. 



In all cases, we use for the ionizing star cluster a 
Salpeter initial mass function with an upper mass limit 
of 120 Mq and a lower mass limit of 0.8 Mq, and we 
assume an instantaneous burst. Sequences of models are 
built with ages spaced by time intervals of 0.5 Myr starting 
from 0.01 Myr. Unless expressed explicitly, the metallicity 
of the nebula is the same as that of the ionizing cluster, 
and the ionized gas contains no dust. 

Our aim is to reproduce the observational diagrams as 
well as possible, with the simplest physically reasonable 
models. 
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The first step is to model the H ii regions not as static 
gaseous spheres or thin bubbles with constant radius as 
was done in all the previous works mentioned in the intro- 
duction, but to explicitly take into account the expected 
time variation of the gas density distribution. As a sim- 
ple representation, we chose the model of supernova/ wind 
driven bubble in an adiabatic phase, which consists of a 
thin spherical shell of density n and radius given by the ex- 
pression R = At^/^ cm, where t is the age in Myr. Such an 
approach has already been used in the past for a detailed 
modeling of several giant H II regions, in which the ob- 
served size was also used as a constraint (Garcia- Vargas 
et al. 1997; Maiz-Apellaniz et al. 1999). In the case of 
constant mechanical energy, the value of A is linked to 
the energy injection E and to the density ng of the am- 
bient medium hy A = 5.1 x 10^'^{E/noY/^ (Weaver et al. 
1977). By trial and error, we take A = 2x 10^° cm which, 
for a typical value of no ~ 1 cm~^ (see e. g. van Zee et 
al. 1998b) corresponds io E ~ % x 10^'' erg s^^. Note 
that, for the stellar population synthesis models we use, 
the mechanical luminosity from stellar winds and super- 
novae during the first 3 Myr of a burst with initial mass 
of 10^ M0 is estimated to be ~ 2 x 10^* erg s"^ at 0.2 Zq 
and 5 x lO^'^ erg s"! at 0.05 Z© (Leitherer et al. 1999). The 
models are built with a density n = 100 cm""^. The inte- 
grated initial mass of the stellar population, M*, is taken 
to be 10^ Mq, which corresponds to a total H/3 luminosity 
of 3.5 x 10"^® erg s~^ at an age of 1 Myr if all the ionizing 
photons are converted into Balmer line photons under case 
B. Note that results from photoionization models (line ra- 
tios and H/3 equivalent width) remain identical to the ones 
presented here for other values of M^, if R is changed ac- 
cordingly (i.e. multiplied by (M^/10''^ M^f-^ so that the 
ionization parameter U = Qn/ (AirR'^nc) is conserved). 

Starting from this simple expanding bubble model, we 
proceed by modifying the model assumptions step by step, 
until the theoretical sequences match all the observed di- 
agrams satisfactorily. Note that, at an age of 0.01 Myr, 
our parametrisation implies a bubble radius of about 4 pc, 
which is the typical size of a stellar cluster. It can easily 
be shown that, for such a situation, our models will over- 
estimate the ionization parameter since they assume that 
all the stars are located at the center of the bubble, while 
the real situation is closer to that of a cluster of Stromgren 
spheres. 

In the forthcoming figures, the models pertaining to a 
given evolutionary sequence are represented by the same 
symbol, and two consecutive models are linked with a 
straight line. Since the models are built at equally spaced 
time intervals, the distribution of observational points in 
a diagram should approximately follow the distribution of 
the symbols. Models having [O III] A4363/II/3 < 2 x 10"^, 
which roughly corresponds to the lower limit of the in- 
tensities of the [O III] A4363 line in our sample, are rep- 
resented by small symbols. Table 1 summarizes the main 
characteristics of the models shown in the figures. 



Table 1. Description of the model sequences 



name 


^* 


^neb 


old stars 


a J 




X-rays 


HI 


0.2 


0.2 


no 


1 




no 


H2 


0.2 


0.2 


yes 


1 




no 


H3 


0.2 


0.2 


yes 


0.5 exp(- 


-t/3) 


no 


T1 

ii 


u.uo 


u.uo 


no 


1 
1 




no 


12 


0.05 


0.05 


yes 


1 




no 


13 


0.05 


0.05 


yes 


0.5 exp(- 


-t/3) 


no 


14 


0.05 


0.05 


yes 


0.5 exp(- 


-t/3) 


yes 


15 


0.05 


0.2 


yes 


0.5 exp(- 


-t/3) 


yes 


L4 


0.02 


0.02 


yes 


0.5 exp(- 


-t/3) 


yes 


L5 


0.02 


0.2 


yes 


0.5 exp(- 


-t/3) 


yes 



4. Results 

4.1. The "high" metallicity bin 

Figure 3 shows, in the same diagrams as Fig. 2, our "high" 
metallicity subsample, which is composed of 105 objects. 
Overplotted with circles is the starting sequence of mod- 
els, with metallicity Z = 0.2 Zq and A = 2 x 10^° cm 
(sequence HI). This sequence of models reproduces the 
classical emission line ratio diagrams [O ill] A5007/H/3 vs. 
[O ii] A3727/H/3 and [O ill] A5007/H/3 vs. [O i] A6300/H/3 
(panels d and e) reasonably well. It also reproduces the 
tendency for [O i] A6300/H/3 to increase as EWQif]) de- 
creases (panel c) due to the fact that, in the expand- 
ing bubble scenario, the ionization parameter strongly 
decreases with time (while in the static models consid- 
ered in the previous studies, it decreases only because of 
the gradual disappearance of the most massive stars). In 
this scenario, there is no difficulty in reaching the highest 
[O i] A6300/II/3 observed. However, this sequence of mod- 
els does not reproduce the slope of the observational points 
in the [O i] A6300/H/3 vs. EW{Ill3) diagram. More impor- 
tantly, it docs a very bad job for the [O ill] A5007/H/? vs. 
EW(Hf3), [O III] A5007/io ii] A3727 vs. EW(Rl3), and 
He I A5876/H/3 vs. EW{Ill3) diagrams (panels a, b and g, 
respectively). As stated in SSLOl, most H ii galaxies are 
likely to contain stars from previous star forming episodes 
which do not contribute to the H/3 flux and reduce the 
value of EW{}i(3) with respect to a pure starburst model. 
Such populations of older stars (hereinafter referred to as 
old populations for simplicity) have been directly inferred 
from deep CCD images (Loose & Thuan 1985; Papaderos 
et al. 1996) and from weak stellar features detected in 
the continuum (e.g., Kong & Cheng 1999; Mas-Hesse & 
Kunth 1999; Raimann et al. 2000; Guseva et al. 2001). 
Statistically, the total flux in the continuum arising from 
old stellar populations should not depend on the age of 
the most recent burst of star formation which provides the 
ionizing stars of the H II regions. For an illustration, we 
show in Fig. 3 as squares the locations of the same model 
sequence as our starting one but with an old population 
whose continuum at H/3 is equal to the stellar continuum 
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from the most recent starburst at the age of 0.01 Myr (se- 
quence H2). Such a small value seems appropriate for our 
objects since the observed spectral energy distributions 
(SED) of young star forming regions in some galaxies from 



the Byurakan sample can be fitted by the theoretical SEDs 
with large fraction of the light from the 3-5 Myr old stel- 
lar populations and minor contribution of the older non 
ionizing populations (e.g. Papaderos et al. 1998; Thuan 
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et al. 1999; Guseva et al. 2001). This brings the values of 
EW(R/3) for the earliest stages in closer agreement with 
the maxinmni observed values. However, this still docs not 
reproduce the slopes defined by the observational points in 
the [O III] A5007/H/3 vs. EW{Rl3) and [O iii] A5007/[O ii] 
A3727 vs. EW(Hi3) diagrams (panels a and b) and the fact 
that He i A5876/H/3 is observed to be constant down to 
EW(R/3) as small as 10. Chosing a larger contribution of 
the underlying population docs not help to explain the 
behaviour of the observational points in panels a, b and 
c. Other effects have been invoked in the literature that 
may reduce EW(Hf3). One of them is the fact that in star- 
burst galaxies, the nebular lines seem to come from dustier 
regions than the stellar continuum (Calzetti et al. 1994, 
Gordon et al. 1997). However, in our sample, the nebular 
extinction is generally very small (see Fig. Ic) so that the 
effect, if existing, will be negligible for the vast majority 
of our objects. Similarly, assuming that a typical aperture 
correction a should be applied to the models will roughly 
translate the curves towards the left, but will not much af- 
fect the slope in these diagrams. What is needed is to ob- 
tain models in which EW{lif3) decreases with time more 
rapidly, without affecting the [O iii] A5007/H/3 or [O iii] 
A5007/[O ii] A3727 ratios. This can be achieved by assum- 
ing that the covering factor / decreases with time, so that 
an increasing fraction of ionizing photons is not processed 
into H/3. Such a scenario is not unlikely, since the true 
geometry may not be spherical but rather axisymmetri- 
cal. Then, during late stages, the matter from the polar 
directions will have been blown out at high velocities, and 
the matter giving rise to the bulk of the emission would 
rather have the form of a ring. This is indeed the geome- 
try suggested from the Ha map of I Zw 18 (Stasihska & 
Schaerer 1999). Another way to obtain a covering factor 
decreasing with time is in the case of gradual clumping 
due to instabilities. To make things simple, we assume 
that time dependence of the covering factor is exponen- 
tial. By trial and error, we find that / = exp(— 1/3), where 
t is the time expressed in Myr, provides a reasonable fit 
to the observed data in panels a, b, c and g, if we main- 
tain the contribution of the old stellar population as in 
the sequence represented with squares and if we assume 
an aperture correction a=0.5. This is likely a conservative 
estimate of the aperture correction in our sample. For the 
closest objects, the Ha emission is more extended that 
the part entering the slit. Therefore, EW{Ill3) as mea- 
sured from spectra is smaller than computed by models. 
Guseva et al. (2000) have determined the aperture correc- 
tion for 33 galaxies from the Byurakan sample and found 
a mean value of 0.33. For distant objects, the slit includes 
the entire H ii region and also an extended old population 
of stars from the underlying galaxy, if present. However, 
as noted in Sect. 2, the fact that there is no correlation 
between redshift and EW{R(3) suggests that, statistically, 
the contribution of the old population to the H/3 contin- 
uum is not very large. The resulting sequence is repre- 
sented by triangles in Fig. 3 (sequence H3). This sequence 
also reproduces the ([O ii] A3727-F[0 iii] A5007)/H/3 vs. 



EW{B.0) diagram. Of course, the inclusion of an old stel- 
lar population and of a decreasing covering factor does not 
change the location of the model points in the line-ratio 
diagnostic diagrams (panels d and e). 

Therefore, we reach the conclusion that the scenario 
of an adiabatically expanding bubble ionized by a coeval 
cluster of stars reproduces both the line ratio diagrams 
and the evolutionary diagrams remarkably well, provided 
that account is taken for an old stellar population and for 
a covering factor decreasing with time. 

The dispersion of the observational points in Fig. 3 is 
perfectly accounted for (within two sigmas) if the value of 
A defining the radius of the bubble lies within ± 0.3 dex 
of the value adopted for the models shown in the figure. 

We note that the He ii A4686/H/3 vs. EW{Rl3) dia- 
gram is however not well reproduced by such a sequence. 
Indeed, He ii A4686 lines are observed at equivalent widths 
between 300 and 50A in our "high" metallicity subsam- 
ple, while the models predict He ii A4686 line intensities 
at the level of 0.01 of H/3 only at EW(iip) around 30A, 
when Wolf-Rayet stars are present. Of course, this value 
of EW{llf3) at which Wolf-Rayet stars appear is depen- 
dent on the underlying stellar population and the cover- 
ing factor, so that some dispersion is allowed. However, 
one does not expect more than 25 - 35 % of cases with 
He II A4686/H/3 of the order of 0.01. Our subsample con- 
tains 10 such objects out of 28 in which the quality of 
the data is sufficient to detect the He ii A4686 line if it 
were present. Of course, one can argue that we are dealing 
with small number statistics and that statistical effects are 
particularly important when considering the He II A4686 
lines which arc produced by only a small fraction of stars 
(Cerviiio et al. 2000). The problem will be clearer when 
considering the next subsamples. 

Also, our models do not reproduce the observed ten- 
dency of [N II] A6584/[0 ii] A3727 to increase as £W(H/3) 
decreases. We will come back to this below. For the mo- 
ment, let us simply note that our former conclusions are 
not affected by this failure. 

4.2. The "intermediate" metallicity bin 

Figure 4 shows the same diagnostic diagrams as Fig. 3 for 
OUT "intermediate" metallicity subsample, which contains 
171 objects. Overplotted with circles, squares and trian- 
gles are the model sequences II, 12, 13. These sequences are 
identical to sequences HI, H2, H3 respectively except for 
the metallicity which is Z = 0.05 Zq, relevant to this sub- 
sample. We immediately see that the models reproduce 
the observational sequence much less satisfactorily than 
in the "high" metallicity bin. Clearly, no solution can be 
found by simply changing the value of A in the expression 
giving the radius of the bubble, the contribution of the 
underlying old population or the law of variation of the 
covering factor. Whatever solution which would increase 
the [O i] A6300/H/3 ratio and bring it in better agreement 
with the observations would simultaneously decrease the 
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[O III] A5007/H/3 ratio and make the situation worse in 
the [O III] A5007/H/3 vs. EWiK/S) diagram. The model 
sequences are significantly displaced to the lower left of 
the observed sequences in the [O ill] A5007/H/3 vs. [O ll] 
A3727/H/3 and [O lii] A5007/H/3 vs. [O i] A6300/H/3 dia- 
grams (panels d and e) and below the ([O ll] A3727-|-[0 III] 
A5007)/H/? zone (panel f). One way to propose a diag- 
nostic is to say that there is not enough heating in the 
models, so that the ratios of collisionally excited to re- 



combination lines are too small. We have attempted to 
remedy this situation by a variety of means. For example, 
we have considered depletion of metals into dust grains, 
as well as heating by photoelectric effect on small grains 
using the prescriptions described in Stasinska & Szczerba 
(2001): this is by far insufficient, given that the dust-to- 
gas mass ratio cannot exceed the maximum available mass 
of elements that can be condensed into grains, and that 
the ionization parameter is very low {U = 10~^ at t = 
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1 Myr, U = 4 X 10~^ at t = 4 Myr). An inhomogeneous 
density distribution does not help either. Obviously, one 
cannot charge the Costar or the Schniutz ct al. (1992) 
model atmospheres used for the spectral energy distribu- 
tion in our models, since they precisely give harder fluxes 
than the more accurate model atmospheres of Pauldrach 
et al. (2001) or Hilher & Miller (1998). 

One clue might come from the He II A4686/H/3 dia- 
gram. Using the models described above, we predict that 
He II A4686 should be measurable (He ii A4686/H/? ~ 
0.01) in less than 20% of objects, (the time resolution of 
our computations actually makes us miss the highest He il 
A4686/H/? value), corresponding to the lifetime of massive 
Wolf-Rayet stars. However, in the Byurakan "intermedi- 
ate" metallicity subsample there are 34 objects out of 62 
with He II A4686/H/3 at the level of 0.01. Therefore, the 
statistical fluctuation argument cannot be invoked here. 
Nuclei of planetary nebulae and hot dwarfs born in pre- 
vious episodes of star formation are not expected to pro- 
vide sufHcient ionizing photons to significantly increase 
the He II A4686/H/3 ratio. Cervifio (2000) and Cervino et 
al. (2002) have invoked X-rays to produce the He ii A4686 
line in giant H II regions. Such X-rays can be produced by 
a variety of sources: individual stars (their contribution 
is in fact considered negligible by these authors), binary 
stars, supernova remnants and hot diffuse gas. In the evo- 
lutionary population synthesis models of Cervifio et al. 
(2002), the X-ray production is estimated as a function of 
an efficiency conversion factor from the mechanical kinetic 
energy injection from stellar winds and supernova explo- 
sions into X-ray radiation. From their study, assuming an 
efficiency conversion factor of 20%, the amount of X-ray 
energy released by a cluster of metallicity 0.05 Zq and 
initial mass of 10^ Mq goes from loss than 10"^* erg s~^ 
to a maximum of 10"^^ erg s~^, which occurs at an age of 
around 3.5 Myr. This is far from sufficient to reproduce the 
observed He ii A4686/H/? ratios, which range between 0.01 
and 0.03 for any value of EW{Ii(3) below SOOA. Assuming 
that the Ho II A4686 ciniission is more concentrated than 
the H/3 one and taking into account that the slit does 
not cover the entire H li region will shift the observational 
points downwards by a factor about 2 (Guseva et al. 2000). 
This is not enough to put into agreement observations and 
model predictions. 

We have construct(xl a sequence of models with the 
same characteristics as the one presented in Fig. 4 but 
with an additional source of X-rays, radiating as brem- 
strahlung at a temperature of 10^ K and having a total 
luminosity of 4 x 10^° erg s~^. This sequence of models 
is shown in Fig. 5 with circles (sequence 14), overplotted 
on the observational data similarly to Fig. 4 (we adopt 
the same contribution of an old population and the same 
law for the covering factor / as previously). From the fig- 
ure, we see that such a luminosity is required to reproduce 
the upper envelope of the observational points. Note that 
these X-rays have little incidence on the other line ratios, 
even on [O l] A6300/H/3 which is enhanced by a factor 2 
at most in the zone of interest, as can be seen by compar- 



ing Figs. 4 and 5. From the range of observed values of 
He II A4686/H/3 we conclude that if these ratios are to be 
explained by photoionization due to the radiation from a 
thermal plasma at around 10^ K, the total required X-ray 
luminosity is 10^° — 4 x 10^'^ erg s~^ for at least half of 
the sources. In any case, we note that photons from such 
an additional X-ray source do not heat the gas strongly 
enough to bring models into agreement with observations 
in the [O III] A5007/H/3 vs. [O ii] A3727/H/3 and [O III] 
A5007/H/3 vs. [O i] A6300/H/3 diagrams (panels d and e). 

A popular way to explain discrepancies between pho- 
toionization models and observed line intensity ratios 
in H II regions (essentially the occurence of high [O i] 
A6300/H/3 or [S ii] A6716, A6731/H/3), is to invoke shocks. 
However, we note that, here, the [O i] A6300 line intensities 
are entirely explained by pure photoionization models of 
reasonable geometry (the same can be said for [S ii] A6716, 
A6731 fines). As a matter of fact, the large [O i] A6300/H/3 
or [S ii] A6716, A6731/H/3 ratios in shocks are mainly due 
to photoionization from X-rays produced behind the shock 
front. Shocks can also heat the gas, and arc sometimes in- 
voked to explain the [O ill] A4363/[0 iii] A5007 ratios 
observed to be higher than obtained in tailored photoion- 
ization models of H ii regions. They could then perhaps 
also explain the [O iii] A5007/H/3 vs. [O ii] A3727/H/? 
diagram, by boosting the collisionally excited lines. The 
mechanical energy from stellar winds and supernovae from 
the latest burst of star formation is actually insufficient to 
affect appreciably the luminosities of the [O III] A5007 and 
[O Ii] A3727 lines, but one could argue that kinetic energy 
is available from previous bursts or from cloud-cloud colli- 
sions. However, it is important to note that what appears 
as a heating problem may actually have a completely dif- 
ferent origin. For example, Stasihska (1999) and Moy et 
al. (2001) have shown that composite models including a 
diffuse component can displace to the left the location of 
model sequences in [O ill] A5007/H/3 vs. [O ii] A3727/H/3 
or [O III] A5007/H/3 vs. [O i] A6300/H/j diagrams with- 
out the need of any additional heating. Chc^niical inhomo- 
genoities can be another appealing explanation. Take for 
example the same sequence of models as shown in Fig. 5 by 
circles, but assume that the metallicity of the gas is Zneb 
= 0.2 Z0(sequence 15), not 0.05 Zq. Then, combine these 
two sequences of models, so as to reproduce an H II region 
which would consist of parcels of gas with composition 
identical to the initial composition of the ionizing stars, 
and of blobs of enriched matter. Such a sequence of mod- 
els is shown in squares in Fig. 5 (here, we have assumed 
that the relative weight of the high metallicity component 
with respect to the low metallicity one, in terms of H/3 
luminosity, is given by the expression exp(— 0.5/t) where 
t is the time in Myr). The [O III] A5007/H/3 vs. [O ll] 
A3727/H/3 and [O iii] A5007/H/3 vs. [O i] A6300/H/3 dia- 
grams (panels d and e) are now quite well reproduced. This 
is also the case for the ([O ii] A3727 + [O iii] A5007)/H/3 
diagram vs. EW(HP) diagram (panel f). This is encour- 
aging, since ([O ii] A3727 + [O iii] A5007)/H/3 is a first 
order indicator of the total energy released in the colli- 
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sionally excited lines, and the basis of abundance deter- 
minations using strong lines only (see Pilyugin (2000) and 
references therein). Of course, we do not mean to say that 
this precise combined model corresponds to reality. There 
are many combinations that can provide an acceptable 
solution from the point of view of line ratios. One must 
be aware that, in such a hypothesis, the abundances de- 
rived by classical methods only represent some average 



abundance, and there is a bias which depends on the ex- 
act abundance pattern (and perhaps also on the density 
pattern). A future work would be to pinpoint which so- 
lutions could be compatible with the theory of produc- 
tion and mixing of elements in dwarf galaxies. A more 
straightforward evidence for self-enrichment of the gas in 
H II galaxies comes from the [N ii] A6584/[0 ii] A3727 
vs. EW{li(3) diagram (panel i). SSLOl already invoked 



12 



G. Stasiriska and Y. Izotov: Modeling the emission line sequence of H ii galaxies 



the possibility of self-enrichment in nitrogen from the be- 
haviour of H II galaxies in such a diagram, irrespective of 
their metallicity. However, they also made the point that 
a similar behaviour could arise from a contribution of old 
stellar populations increasing with metallicity. In all the 
figures of the present paper (except Fig. 2) the objects are 
grouped by metalhcity, and the trend of [N il] A6584/[0 ii] 
A3727 increasing with decreasing EW{H.fi) remains. It is 
even more conspicuous in the intermediate metallicity bin 
than in the high metallicity bin (which, remember, con- 
cerns only objects with 0/H < 2 x 10~^, i.e. much smaller 
than the solar value). 

4.3. The "low" metallicity bin 

Only a small proportion of objects from our sample have 
0/H determined to be smaller than 3 x 10~^: 15 in to- 
tal. Indeed, it is a known fact (e.g., Terlevich et al. 1991; 
Masegosa et al. 1994; Kunth & Ostlin 2000) that H ii 
galaxies with metallicities smaller than 1 /25 solar are very 
rare. Therefore, the trends in observational diagrams arc 
less clearly seen than in the case of the higher metallic- 
ity bins. Nevertheless, by comparing the location of the 
observational points in Fig. 6 with those of Fig. 5, we 
note that the same trends are present. Superimposed on 
the data points, we represent with circles the model se- 
quence L4, which is identical to the model sequence 14 
except for the metallicity which is Zneb = Z^, = 0.02 Zq 
instead of 0.05 Z©). We see that this sequence is a good 
first approximation to the observed trends (as was the 
case with sequence 14 for the intermediate metallicity bin) . 
But, similarly to the intermediate metallicity bin, the pre- 
dicted intensities of coUisional lines with respect to H/3 
are lower than observed. The problems seems even exac- 
erbated: compare the position of the circles with respect 
to the data points in Figs. 6 and 5. This can be qualita- 
tively understood if self-pollution is indeed the explana- 
tion: the effects of self-pollution are expected to be larger 
for lower initial metallicities. For an illustration, we plot 
with squares a composite model sequence similar to the 
one plotted in Fig. 5, except that Z^eb is 0.02 Zq and 
0.2 Zq in the low and high metallicity components re- 
spectively and that is 0.02 Z©. 

5. Discussion 

Using a sample of H ii galaxies of unprecedented size 
(about 400 objects) in which the oxygen abundances have 
been obtained using the temperature derived from the 
[O III] A4363/5007 line ratio, we have confirmed the im- 
pressive trends and extremely strong correlations found in 
diagrams relating line ratios and H/3 equivalent width that 
we already found in previous studies (SL96 and SSLOl). 
The He II A4686 line is present at a level of 1 - 3 % of 
H/3 in a large number of H II galaxies, roughly half of the 
galaxies where the signal-to-noise was sufficient to detect 
and measure this line (which often appears on top of a 
Wolf-Rayet bump in the spectra). This line is present in 



the entire range of H/3 equivalent widths, and its intensi- 
ties show no correlations neither with EWilipi) nor with 
metalhcity. The increase of [N Ii] A6584/[0 ii] A3727 with 
decreasing EW{Yi(i), already noted by SSLOl, is clearly 
seen (although more dispersed that some other trends). 

So far, analysis of emission line trends in H ii galaxies 
made use of grids of photoionization models, but without 
attempting to propose a sequence of evolutionary models 
reproducing the observed trends. 

Taking advantage of the large size of our sample and 
of the fact that the metallicities were determined in a di- 
rect and independent way, we divided our sample in three 
metallicity bins. We considered each bin separetely, try- 
ing to find under what conditions a sequence of photoion- 
ization models, taking into account the time evolution of 
the ionizing cluster, reproduces all the observational dia- 
grams adequately. In such an exercise, we used only very 
simple prescriptions, justified if possible by arguments re- 
lated with what is known of the physics of these objects. 

We have found that the simple model of an adiabatic 
expanding bubble gives rise to an increase of the [O i] 
A6300/H/3 ratio with time, reaching the highest values ob- 
served in about 5 Myr. The presence of an older, non ion- 
izing stellar population that contributes to the continuum 
at H/3 has been attested by direct studies of the stellar 
features in the continuum of H II galaxies. The charac- 
teristics of such an old population should of course be 
independent of the age of the most recent burst of star 
formation. Since the contribution of this old population 
to the stellar continuum at H/3 is proportionally larger for 
younger starbursts, such an underlying population is not 
sufficient to reproduce quantitatively the observed trends 
in the [O III] A5007/H/3 vs. EW{Y{f3), [O ill] A5007/[O ii] 
A3727 vs. EW(}il3), and He i A5876/H/3 vs. EW{}l(i) di- 
agrams. The observational data require that the covering 
factor / of the ionizing source by the emitting gas de- 
creases with time. A reasonable fit to the observations is 
provided assuming an aperture correction of 0.5 to the 
models and an exponentially decreasing covering factor, 
with an e-folding time of 3 Myr, associated with an old 
population contributing at a level equal to the H/3 contin- 
uum produced by the young stars at zero age. That the 
covering factor of giant H ii regions is smaller than unity 
and depends on the evolutionary stage of the region has 
already been suggested by Castellanos et al. (2002) from 
the study of a limited number of objects in spiral galaxies. 
What we find here is that there seems to be an universal 
law for the temporal variation of the covering factor. The 
physical reason underlying such a law could be linked to 
departure from spherical symmetry and blowout or to shell 
fragmentation (see Tenorio-Tagle et al. 1999). 

We find that the He ii A4686 nebular line emission in 
H II galaxies occurs too frequently and in too wide a range 
of EW{i{j3) to be attributable to either the hard radiation 
field from Wolf-Rayet stars, as was suggested by Schaerer 
(1996, 1998), or the X-rays produced by current evolution- 
ary synthesis models of young starbursts (Cervino et al. 
2002). Assuming that the He ii A4686 line is due to pho- 
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Fig. 6. The data points correspond to the "low" metallicity bin, The data points have 0/H < 3 x 10~^. The model 
sequence represented by circles corresponds to sequence L4: metallicity Z = 0.02 Zq, expanding bubble with varying 
covering factor and old underlying population included, and additional ionization by X-rays. The model sequence 
represented by squares corresponds to the composite sequence L5, i.e. same as L4 but including the effect of self- 
enrichment (for complete descriptions of the model sequences, see text and Table 1). 



toionization by a hot plasma at a temperature of 10^ K, a 
total X-ray luminosity of 10^" — 4 x 10^" erg is required 
for at least half of the sources. In some of the objects, 
the observed He li A4686 nebular emission could actu- 
ally be due to harder X-rays produced either by massive 
binaries or by supernova remnants from previous genera- 
tions of stars in the age range of 10 - 50 Myr (Van Bever 



& Vanbeveren 2000). Upcoming X-ray spectroscopy and 
imaging will help clariiy the matter in the near future. 

Our evolutionary sequences of models fail to repro- 
duce perfectly the [O iii] A5007/H/3 vs. [O ii] A3727/H/? 
and [O III] A5007/H/3 vs. [O i] A6300/H/3 diagrams for 
the intermediate and low metallicity bins, if the chemi- 
cal composition is homogeneous and identical to the one 
of the gas which gave birth to the young stars. The col- 
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lisionally excited lines are slightly - but significantly - 
too weak with respect to the observations. This problem 
had been noted before (SSLOl), but here it is more visi- 
ble due to the fact that we compare models with objects 
of similar metallicity. One way to understand the discrep- 
ancy is to postulate the existence of an additional heating 
mechanism. It is unlikely that this heating agent could be 
dust, or shocks produced by the energy released by winds 
and supernovae from the latest burst of star formation. 
Planetary nebulae and hot white dwarfs from earlier gen- 
erations of stars would also be far from sufficient. The hy- 
pothesis of shock heating due either to stellar winds from 
previous generations of stars, or to cloud-cloud collisions 
needs to be investigated. However, another option to ex- 
plain these emission line diagnostic diagrams is to invoke 
chemical inhomogeneities. Chemical inhomogeneities are 
expected in regions experiencing mass loss and supernova 
explosions from massive stars. The question is in what 
form is the newly synthesized matter, and how much of 
it has escaped the nebula. The problem is far from be- 
ing settled at present (Roy & Kimth 1995; Tenorio-Tagle 
1996; Silich & Tenorio-Tagle 1998; Mac Low & Ferrara 
1999; Ferrara & Tolstoy 2000). So far direct evidence for 
self-enrichment has been scarce (see Kobulnicky (1999) 
for a review), but the fact that [N ii] A6584/[0 ii] A3727 
increases as EWCRP) decreases, especially in the interme- 
diate and low metallicity bins of our sample, argue for the 
existence of self-enrichment in nitrogen on timescales of a 
few Myr. Of course, this in itself does not argue for any 
self-enrichment in oxygen, since the nitrogen and oxygen 
donors are different (see Matteucci & D'Ercole (1999) for 
a review). Given that the fate of the elements is still not 
well understood, we have felt it instructive to investigate 
the effect of inhomogencous metallicity on the emission 
line diagnostic diagrams. We found that for the interme- 
diate metallicity bin, a combination of two models with 
abundances Zneb ~ 0.2 Zq and 0.05 Zq leads to good 
agreement with the observations, but other solutions are 
likely possible as well. For the low metallicity bin our mod- 
els indicate that a qualitatively similar composite model 
would be acceptable. The scatter of observational points 
in that bin is larger than in the other bins. This can be 
attributed to the small number of points in that bin and 
also to the fact that the effects of self-enrichment become 
more drastic at low metallicity. It is only for our high 
metallicity bin that there is a priori no need to invoke self- 
enrichment in oxygen, since the diagrams are satisfactorily 
reproduced for models with homogeneous chemical com- 
position. However, chemical inhomogeneities, if present, 
are expected to have a smaller impact, due to the fact 
that [O III] A5007/H/3 and [O ii] A3727/H/3 are not so de- 
pendent on metallicity for Zneb between 0.2 - 1 Zq (see 
e.g. Stasihska 2002). 

Finally, it is worth of noting that our prescriptions 
for the time dependence of the covering factor and of 
the degree of inhonomogeneity were tailored to rc;prodiice 
the emission line sequence of H ii galaxies. On the other 
hand, the time dependence of the radius of the expand- 



ing bubble is exactly the one predicted by the theory of 
supernova- and wind-blown bubbles. The size of the bub- 
ble required to reproduce the observational diagrams cor- 
responds rather well with the prediction using the energy 
available from stellar winds and supernovae computed 
in population synthesis models. It is most encouraging 
that emission line diagnostics of H ii galaxies confirm, on 
a statistical basis, the predictions of theories on stellar 
evolution, stellar populations and dynamical interaction 
with the interstellar medium. However, our interpretation 
should be compared to detailed modeling of selected ob- 
jects in the entire sequence of EW(H(3). In particular, it 
would be important to compare the ages derived directly 
from the analysis of the stellar UV light with those implied 
by our interpretation of the emission line diagrams, and to 
use Wolf-Rayet stars as independent clocks. So far, only a 
handful of galaxies from our sample have been studied in 
such a way (Mas-Hesse & Kunth 1999), and the compari- 
son is not conclusive. 
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